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Abstract 

In this paper we discuss a minor modification of a previous SU(5) x A 5 flavour model 
which exhibits at leading order golden ratio mixing and sum rules for the heavy and the 
light neutrino masses. Although this model could predict all mixing angles well it fails 
in generating a sufficient large baryon asymmetry via the leptogenesis mechanism. We 
repair this deficit here, discuss model building aspects and give analytical estimates for 
the generated baryon asymmetry before we perform a numerical parameter scan. Our 
setup has only a few parameters in the lepton sector. This leads to specific constraints 
and correlations between the neutrino observables. For instance, we find that in the model 
considered only the neutrino mass spectrum with normal mass ordering and values of the 
lightest neutrino mass in the interval 10 — 18 meV are compatible with the current data 
on the neutrino oscillation parameters. With the introduction of only one NLO operator, 
the model can accommodate successfully simultaneously even at Icr level the current data 
on neutrino masses, on neutrino mixing and the observed value of the baryon asymmetry. 
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1 Introduction 


The theoretical explanation for the observed neutrino oscillations and neutrino masses re¬ 
quires physics beyond the Standard Model. Furthermore the presence of Dark Matter and 
the observed baryon asymmetry of the Universe (BAU) support the need for a more funda¬ 
mental theory. In the present article we will establish a connection between two of the above 
mentioned observations and investigate the Baryogenesis through leptogenesis scenario [1] 
in an SU(5) x A 5 flavour model. The model we are going to discuss here is the first GUT 
A 5 golden ratio flavour model with successful leptogenesis to our knowledge. This recently 
proposed model [ 2 ] has the feature that 6*12 is connected to the golden ratio (jig = via 

612 = tan“^ (^)’ golden ratio (GR) type A models in [3] the reactor angle is 

predicted to be vanishing at leading order and the atmospheric angle to be maximal. Hence, 
the neutrino mixing matrix Ugr has the form 


Ugr 


fPg 

V U5 

\l 4>g'P> 

0 1 

/ 1 
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/ 24 , 

V 2V5 

U2 

1 i 



24>gV5 

V 2 U 5 

V2 ) 


Pn 


( 1 . 1 ) 


which is given in the convention of the Particle Data Group [4] with the diagonal matrix 
Pq = Diag(exp(—exp(—^), 1) containing the Majorana phases. Since the experimental 
values for the angles, cf. Tab. 1, strongly disfavour ^13 to be vanishing the leading order 
mixing angles have to be corrected to realistic values. In [2] we followed the approach based 
on Grand Unification where the neutrino mixing angles receive corrections from the charged 
lepton sector. Namely this model features SU(5) unification. Thereby we could explore the 
SU(5) relation 6*13 Ri dclV2 from the non-standard Yukawa-coupling relations yr/yb = —1.5 
and y^/ys = 6 [6,7], and for the double ratio {yg/ys){yd/ye) = 12 which are all in perfect 
agreement with experimental data. 

In addition to the corrections from the charged lepton sector renormalisation group run¬ 
ning effects (RGE) have to be taken into account. Due to a neutrino mass sum rule in both 
hierarchies only a certain mass range is allowed. For the inverted ordering this implies large 
RGE effects for 612 which rule out this ordering. 

Since the light neutrino masses are generated via the type-I-seesaw mechanism in [ 2 ] the 
Baryogenesis through leptogenesis mechanism can be easily implemented. In this mechanism 
the dynamically generated lepton asymmetry is converted into a baryon asymmetry due to 
sphaleron interactions. Thermal leptogenesis can take place when the heavy RH Majorana 
neutrinos (and their SUSY partners the sneutrinos) decay out-of-equilibrium in a GP and 
lepton-number violating way. Flavour effects [8-10] (see also, e.g., [11-14]) can play an im¬ 
portant role in thermal leptogenesis. We set the scale at which leptogenesis takes place to be 
the see-saw scale Ms = 10^^ GeV. In the model considered we have also tan/3 = 30 [2], and 
thus, the scale Ms falls in the interval 10^(1 -|- tan^/3) GeV< Ms < 10^^(1 -|- tan^/3) GeV. 
For values of Ms in this interval [15] the baryon asymmetry is produced in the two-flavour 
leptogenesis regime and we perform the analysis of baryon asymmetry generation in this 
regime. 

We will see that the original model cannot accommodate for the observed value of the 
baryon asymmetry due to the structure of the neutrino Yukawa matrix. In fact, there would 
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Parameter 


best-fit (Tier) 

3cj range 

9i2 in ° 


33 . 4810:^8 

31.29 ^ 35.91 

6*13 in ° 


8.5010;20 

7.85 ^ 9.10 

6»23 in ° 


42.3116 

38.2 ^ 53.3 

(5 in ° 


306lf^ 

0 ^ 360 

in 10 “® 

eV 2 

'-^^-0.17 

7.02 ^ 8.09 

in 10 “^ 

eV 2 

^•^^'-0.047 

2.317^ 2.607 


Table 1 : The best-fit values and the 3cr ranges for the parameters in the normal ordering 
taken from [5]. 


be no baryon asymmetry generated via the leptogenesis mechanism. In order to generate a 
non-zero asymmetry we will introduce only one additional operator in the neutrino sector 
which corrects the neutrino Yukawa matrix and subsequently affects the phenomenology of 
the model. 

The paper is organized as follows: Section 2 is a short overview of the model building 
aspects including the NLO operator. In section 3 we discuss the analytical results for the 
phenomenology of the model. There we also describe the relevant formulas for leptogenesis. 
In section 4 we show the results of a numerical parameter scan. We discuss the predictions for 
the mixing parameters including the phases as well as for the sum of the neutrino masses, the 
observable in neutrinoless double beta-decay, the kinematic mass and for the generated 
baryon asymmetry. In section 5 we summarise and conclude. 

2 Model building aspects 

The model we are going to discuss is based on the SU(5) x A 5 model proposed in [2]. We only 
had to extend it minimally to accommodate successful leptogenesis. The modification we 
are going to introduce has further implications for the phenomenology. We focus first on the 
related model building aspects. We briefly revise the leading order (LO) superpotential for the 
neutrino sector, which is identical to the original model before we introduce the corrections. 
They are induced by an additional operator in the superpotential which yields next-to-leading 
order (NLO) corrections to the Yukawa couplings while the right-handed neutrino Majorana 
mass matrix remains unaffected. This single higher order operator will generate a sufficiently 
large baryon asymmetry to be in agreement with the experimental observations. 

2.1 The neutrino sector at LO 

We briefly summarise next the relevant parts of the original SU( 5 )xA 5 model from [2]. We 
are not going to discuss the flavon vacuum alignment which does not change at all and is 
given in the original paper. 

The matter content of our model is organised in ten-dimensional representations of SU(5), 
Ti with i = 1, 2, 3, five-dimensional representations F, and one-dimensional representations 
N which transform as one-, three- and three-dimensional representations of A 5 respectively. 
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Table 2: Charges under Z„ and SU(5) and A 5 representations of all the fields appearing in the 
neutrino sector of the model. Note that the only new fields compared to [2] are the messenger 
fields 4'i, ^ 1 , ' 1>2 and ^ 2 - 


see also Tab. 2. 

Additionally, we had introduced in the original model the following flavons which will 
appear in the neutrino sector. There are two flavons which transform as one-dimensional 
representations under A 5 

{O 2 ) = ve^ , (ei) = t’ei , (2.1) 

one flavon in a three-dimensional representation 

( 0 , 0 , 1 ) , ( 2 . 2 ) 

and one flavon in a five-dimensional representation 


(w) 



{V2 + V3),V3,V2,V2,V3 


( 2 . 3 ) 


Their charges under the shaping symmetries are given in Tab. 2. Note especially that no new 
flavon appeared. 

We are not going to discuss here the quark sector, it was analysed in [2]. In what concerns 
the charged lepton sector, we only note that for the matrix of charged lepton Yukawa couplings 
we find 


Te = 


— 1/2021 0 
6022 6032 

0 — 3/2033 


0 

6012 

0 


3 


(2.4) 



The order one coefficients in front of the parameters aij are SU(5) Clebsch-Gordan coefficients 
which imply that 

(2,5) 




where 9c is the Cabibbo angle. This is only possible due to the non-standard Clebsch-Gordan 
coeffcients [6] as it was realised in a series of papers [7]. 

The flavon oj is responsible for the GR structure of the Majorana mass matrix which can 
be seen in the LO superpotential for the neutrino sector which reads 


= ViFNH^ + y'^NNoj . 


The right-handed neutrino mass matrix then reads 

-\/3u2 


(2\II{v2 + Vs) 


3' 

-•v/3u2 


Mrr = 

y -\/3U2 

and the neutrino Yukawa couplings are 




^^3 


3(^2 + t’ 3 ) 


-\/3U2 

1 {V2 +^^ 3 ) 

\/6U3 


\ 


yLO _ n 

— Vl 



( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


which are diagonalised by the golden ratio mixing matrix 17 gr from eq. (1.1). Note that we 
are using the right-left convention for the Yukawa matrices, which means that the first index 
of the matrix corresponds to the SU(2)i singlet. 

This is the structure of the original model which cannot accommodate for the observed 
value of the baryon asymmetry of the universe, as we will see later on. In order to generate 
a non-zero asymmetry we follow the approach as described in [16] and introduce an addi¬ 
tional operator which perturbs the original flavour structure of the model. Note nevertheless, 
that compared to [16] we do not introduce an additional flavon and no additional shaping 
symmetries. Here it is sufficient to extend minimally the messenger content of the model. 


2.2 The neutrino sector at NLO 

In this section we discuss the NLO superpotential of the neutrino sector. In order to accom¬ 
modate leptogenesis which can generate the observed value of the baryon asymmetry, we will 
introduce a correction to the neutrino Yukawa matrix governed by the operator FN(f)s 62 Hse\. 
This operator was absent in the original model but can be added by introducing only two 
new pairs of messenger fields Ti, Ti, T 2 and ^ 2 - 

The renormalisable superpotential of the NLO neutrino sector reads 

yyren, NLO ^ + M^,^-^2^2 + TTiTg + Y’Li^2 + ^2775ri + eieif 3 + 02</>3f 1 , 

(2.9) 

where we have omitted the coupling constants to increase clarity and only write down the 
operators which are new. We assume the messenger masses to be larger than the GUT scale 
and to be related to the messenger scale A by 0(1) coefficients. The charges of the new 
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Figure 1: 


The supergraphs for the neutrino sector including LO and NLO operators. 


messenger fields under the symmetries as well as their SU(5) and A 5 representations are 
shown in Tab. 2. The corresponding supergraphs for the leading order operator and the 
next-to-leading order operator for the neutrino Yukawa matrix can be found in Fig. 1. 

No new operators compared to the original model are possible apart from the one we 
discuss now. The only new effective operator is 

, ( 2 . 10 ) 

where we have denoted with brackets the A 5 contractions. This operator gives a correction 
to the neutrino Yukawa matrix 


/O 1 0\ 

6Y, = \y^\cF^ -1 0 0 , ( 2 . 11 ) 

Vo 00/ 

where 0 < c <C 1. This correction disturbs the golden ratio mixing pattern already in the 
neutrino sector by itself and subsequently the phenomenology of the original model, especially 
the prediction for leptogenesis is modified. 

3 Phenomenology: Analytical Results 

In this section we will discuss the phenomenological implications of introducing 51^. Because 
c is small, in many cases the results are similar to those obtained in the original model. 
However, as we will see, in some cases when the leading order result was relatively small, a 
correction of order c can have a sizeable impact. 
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3.1 Masses and Mixing Angles 

The original model was very predictive due to the built-in sum rules. And indeed one mass 
sum rule remains valid. Since Mrr is not corrected the sum rule for the right-handed neutrino 
masses 

Ml + M 2 = Ms, (3.1) 

is still correct. Note that, the masses are taken here to be complex. 

The situation for the light neutrino masses is somewhat different: they get corrections of 
the order of <C 1. However, since these corrections are small, the sum rule 

ei«i 1 

-r I—r ~ 1 — 

mil |m2| |m3 

is still a good approximation. And hence our estimate for the ranges of the neutrino masses 
from [ 2 ] 


= 0 {c^ 


(3.2) 


0.011 eV < mi for NO, (3.3) 

0.028 eV < m 3 < 0.454 eV for 10. (3.4) 

is still reasonable. For all three PMNS mixing angles we find corrections of order c. For 9is 
and 023 the expressions are somewhat lengthy and not insightful, but as an example we find 
as correction for 0 i 2 to first order in c 


CCOS7 ^/ 2 c^{M^M2) sin7 

^71 |M2| - |M|| 


(3.5) 


The corrections have immediate consequences for the phenomenology. 

The first thing one might wonder, is if the inverted mass ordering is still excluded like in 
the original model. We begin our discussion with the sum rule [17]^ 


sin^ 012 


cos^ 05^2 + 


sin 2012 sin 0 i 3 cos 5 — tan 623 cos 20(^2 
tan 023(1 - cot^ 023 sin^ 013 ) 


(3.6) 


where cos^ 612 = and cos 20(^2 = 2(l)g/\/5 — 1 . 

For c = 0 we can evaluate this sum rule and find 0 i 2 ^ 23° (compared to 0 i 2 ^ 24° from 
the sum rule in [18]). Using the lower bound on the mass scale in the inverted ordering case, 
cf. eq. (3.4), we can estimate that the RGE evolved value of 0 i 2 at the seesaw scale has to be 
smaller than about 5.7°. Attributing this difference completely to the correction 6612 with 
7 = 0 or TT we find c > 0.43. This is a crude estimate because the other mixing angles are 
affected as well modifying the above sum rule in a non-trivial way. Nevertheless, from here 
we would still expect c to be of order 0.1 to safe the inverted mass ordering. On the other 
hand such a large value for c is not plausible from a model building point of view because it 
is associated to a highly suppressed operator making values of c = 10 “^ to 10 “^ plausible. 
We will come back to this point later when we discuss the numerical results. 

^In the original model we had used another sum rule from [18] which can be derived from this sum rule by 
expanding in O 13 . But since ^13 is not very small we want to use now the improved sum rule. 
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For future convenience, we introduce the following definitions for the parameters which 
we will use as well in our numerical scan 

Mi = ^(X + y) = ^|X||l + />e'^|e'^S 0i = arg(X + y) , (3.7) 

M 2 = ^{X-Y) = ^\X\\1-p e'^|e'^^ 02 = arg(X - y) , (3.8) 

Ms = = y||X| 03 = arg(y) , (3.9) 


where 


X = (Ms + V2)y2 , 
y = 3V5v2y2 , 

Y 

X ’ 

0 = arg(y) — arg(X) . 


(3.10) 

(3.11) 

(3.12) 

(3.13) 


In this way, we express the absolute value of the three heavy neutrino masses in terms 
of three real parameters, i.e., |X|, p and 0. \X\ sets the scale of our interest, p reflects the 
detailed structure of the heavy neutrino mass spectrum. 0 is connected to p via the ratio of 
two mass squared differences 


Am^i 16/9 cos 0 

Am|i (/O^ — 2/9COS 0 + 1)(/92 + 2pcos 0 — 3) 


(3.14) 


Notice that we neglect here for the moment RG effects on the masses and corrections of order 
(? to the neutrino masses which will turn out to be well justified in the numerical analysis. 

One of the Majorana phases which we choose to be 0i can be set to zero by applying a 
redehnition of the heavy Majorana fields. The remaining two phases 02 and 03 can as well 
be expressed in terms of p and 0 using the complex mass sum rule Mi + M 2 = M 3 


IM3P - |Mip - IM2P _ l-p^ 

21 M 1 IIM 2 I cos 2 (f> + p^ 

|Mip - IM 2 P + IM 3 P _ 1 + /9COS0 

21 M 1 IIM 3 I “ cos 0 + /9^ 


(3.15) 

(3.16) 


Notice that only normal ordering is viable in this model, and the Yukawa couplings are 
degenerate in LO so that we have \Ms\ < IM 2 I < |Mi|. Thus cos0 is positive, cos 02 is 
negative, cos 03 is positive, which gives us a first constraint on p which we will comment on 
later. Notice also that the sign of sin 02 and sin 03 is not hxed. We plot the dependence of 
the phases on p in Fig. 2. The Majorana phases ai, q ;2 and 02 , 03 are related via 


= -03 and a 2 = 02 - 03 • (3-17) 

up to order (?. 

We comment a little on the phases in the mass matrices. The heavy neutrino mass matrix 
is diagonalised as 

C/;^rMrr[/gr = Dn = Diag(Mi e'^S M 2 M 3 e'^^) ^ Mr 2,3 > 0 . (3.18) 


7 




4^ ^ ^ 5 ^ ^ g ^ ^ 7 ^ ^ 8 ^ ^ ^ 


P 



Figure 2: The dependence of (j) and the phases of the heavy Majorana neutrinos on p according 
to eqs. (3.14, 3.15, 3.16). The unit for </), 02,</>3 is rad. We use the best fit values from the 
global fit [5] for the two squared mass differences as input here. 
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We eliminate the common phase by setting cpi = 0 and attribute the phase factors to a phase 
matrix P = Diag(l, Thus we have 

P-1C/;^rMrr[/grP-i = Diag(Mi,M 2 ,M 3 ) , (3.19) 

which means Ug^P~^ diagonalises the heavy neutrino mass matrix to real and non-negative 
eigenvalues. Applying the seesaw mechanism, we have 

mhL = = U*Bia,g{mi,m2,m3)Ul . (3.20) 

Notice that from Uq^^iti-oU gr = yfwDiag(l, 1, —1) + 0(c) we get 

U, = iUGRP +0(c). (3.21) 

If we choose y" to be real and positive, the only complexity comes from the heavy neutrino 

mass matrix and the phase 7 . By ascribing the phases to the P matrix, the rrii, i = 1 , 2,3 

in eq. (3.20) are real and non-negative up to corrections of 0{c). From now on, we use the 
symbol Mi and rrii to label the real and non-negative masses. 


3.2 Leptogenesis 

In this section we discuss analytical estimates for the generated baryon asymmetry including 
all relevant parameters and formulas. To discuss leptogenesis in this model, we hrst set 
our scale of interest, the see-saw scale, to be Ms — 10^^ GeV. Taking into consideration 
tan/3 = 30, we have 10^(1 -|- tan^ /3) GeV< Ms < 10^^(1 -|- tan^ /3) GeV, which as was shown 
in [15], corresponds to the “two-flavoured leptogenesis” regime [8-10], i.e., the regime where 
the processes mediated by the r Yukawa couplings enter into equilibrium. Later on in our 
numerical scan we will find that in order to generate realistic neutrino masses the parameter 
p has to satisfy the inequality p > 5.8. This in turn implies that the leptogenesis regime 
in the model we are considering cannot be resonant. Indeed, as can be shown, for p > 3.7 

we have Mi — M 2 ;§> Fi = ^ the smallest heavy Majorana neutrino mass splitting 

(Ml — M 2 ) being by at least two orders of magnitude larger than Fi. Thus, the condition of 
resonant leptogenesis [19] (Mi — M 2 ) ~ Fi, is not satished in the model under discussion. 

The CP-asymmetry generated in the lepton charge L; by neutrino and sneutrino decays, 
I = e,pL, T, is [13]; 



1 1 


j;im[(W),Kh;)*i(h;^J)ii]/(((|) 


+ 


1 1 

^{Y,Yj)ii 


Y,M(Yu)ji(Y.)tiiY,Yj),,] 


m. 


2 2 
m, — mf 

j * 


(3.22) 


where 

/(*) =-^ (^ 2 ^+^1 +^)) • ( 3 - 23 ) 

The second term in eq. (3.22) corresponds to the self-energy diagram with an inverted 
fermion line in the loop. It would vanish when we sum over a, and we would end up with the 
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same formula as in the one flavour case 


e- = Ve' = — V f(rn^) 

* Y (y^yh^ 


(3.24) 


In the basis where the charged lepton and the right handed neutrino mass matrices are 
diagonal, we have 

Y, = {UGRP-^)Hyu + SY,)Ue 

( e''’l2+V3+v^cfo') e-'^l2 ( c= e‘'^l2 - Vs+v^sfo) , \ 


= - 


^ ((5+v^)e‘'^i2sf2-2vdI)cf2) e^"‘'^l2 (^(5+V5)cf2e‘'*i2+2vd0^ 

2J5{5+V5) 2 53/4^/1+71 




^^ 1 ( 25 ^ 2 + 03 ) 56 ^ 


V 2 

(cf2-V3+V5e-^l2 5f2) 


'^12 (73+\/5cf2e‘'^l2 +^ 12 ) 


n ei7 (2v^sf2e‘"l2+75cf2+cf2) e^--^i2 (-2y2cf2e-"l2+75^f2+sf2) 


(3.25) 


where we use Ue — Uu, and the abbreviations sin 6*^2 = '§12 and cos 0^2 = 02 - 

Here and in the following we have used the freedom to redefine 1/ by a global phase to 
make y” > 0 so that we find 


F,Fj = PU^j,{Y,Yj + YJY^ + 6Y,Y^)UgrP-^ 

( 

= (y^fl + ciy^f \ -iv^sin7e‘^2/2 


(3.26) 


iv^sin 7 e-‘'/ 2/2 


■ , 3 cos7e~‘('^3-'/2)/2 

V5+y5 


y-yi+/gCOS7e‘'^3/3 _ ^ ^2^ cos7e‘('/3-</2)/2 0 J 

which we have expanded up to 0 {c). 

We give next the expressions for the CP-violating asymmetries in the I lepton charge Li, 
generated in the decays of the heavy Majorana neutrinos A^i, N 2 and N 3 , as calculated from 
eq. (3.22); 


c{y^f 1 


e[ = —7-7= sin7cos(/>2/(^) -siny^-2 

oTT VlO V nT ,2 — rnf 


+ cos 7 sin 


cjy^f 1 
Stt 710 
Cjy^f 1 
Svr 7 T 0 


1-2 "0-^ 
^2 


(3.27) 


2 \ 

-sin 7 cos(^ 2 /(^)+sin 7 ^ - cos 7 sin ((^3 - 02)/(^) ) , (3.28) 

nil '^2 / 


cos 7 (- sin 03/(^) + sin ((/.g - (p 2 )fi^)) 


(3.29) 


We see that to leading order, = 0 and hence leptogenesis was not viable in the original 
model. As 7, in leading order is unitary (except for an overall factor (yf)^), we have e| = 
ef + ef = —ej to leading order. 
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Notice that we compute the CP asymmetry generated by all three heavy (s)neutrino 
decays since the heavy neutrino spectrum is not very hierarchical in our case. 

At the leptogenesis scale and values of y” of interest, the AL = 2 processes are negligi¬ 
ble. They would be important for a different setup with maximal perturbative values of the 
Yukawa coupling of interest (say, for y” = 1) if the leptogenesis scale would be 10^^ GeV (or 
for masses of the heavy Majorana neutrinos of the order of 10^^ GeV). Thus, we can use the 
following analytic approximation for the efficiency factors [10], which accounts for the AL = 1 
interactions and the decoherence effects: 


vimi) 


mu 


.25 X 10-3 eV 


-1 


+ 


0.2 X 10-3 eV 

mi 


- 1 . 16 ' 


-1 


where 


mu = 


vlliY. 


V )il 


Mi 


(3.30) 


(3.31) 


where we introduce another index i, i = 1,2,3, to label the correspondence to the i-th 
heavy (s)neutrino and I = e,pL,T. If we only keep leading order term in rhu, we will have 


+ mu, = - mr- We list the washout mass parameters as follows 


Tni2 — ^26 


1 

5 + V5 Ml 


mir = 


m- 2 r = - 1 + ^ 


1 \ 


\/5 / M 2 






2 M 3 


(3.32) 

(3.33) 

(3.34) 


We do not include the higher order terms 0{c) in the expressions for rhu because they generate 
subleading insignificant corrections. The baryon asymmetry generated by each heavy neutrino 
decay is [15] 


/494 


V761 


mir - v 


V ~ -3 X 10“^e[ ( y 

and the total baryon asymmetry is 

Yb = '^Yi, i = 1,2,3. 


541 

761 


mi2 


(3.35) 


(3.36) 


Notice that we use an incoherent sum over the asymmetry generated by each heavy (s)neutrino. 
This approximation corresponds, in particular, to neglecting the wash-out effects due to the 
lighter heavy Majorana neutrinos iV 2,3 in the asymmetry generated by the heaviest Majorana 
neutrino Ni. Thus, we effectively assume that the indicated wash-out effects cannot reduce 
drastically the asymmetry produced in the Ni decays. Since the masses of ^ 2,3 and Ni in 
the model we are considering differ at most by a factor of 5, we can expect that at least for 
some ranges of values of the masses of Ni and ^ 2,3 the wash-out effects under discussion 
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4 ^ ^ ^ ^ 5 ^ ^ ^ ^ 6 ^ ^ ^ ^ 7 ^ ^ ^ ^ 8 


P 

Figure 3: The single asymmetries Yi,y 2,^3 and the total asymmetry Yb- In the upper plot 
we use c = 0.05 « 6»^, 7 = 2/c7r (fc = 0, ±1, ±2,...), u„ = 246 GeV, = 0.1, |X| = lO^^ QgV 
and in the lower plot we use c = 5.8 • 10“^ « 6^, 7 = tt/ 2 + 2kn, Vu = 246 GeV, y" = 0.1, 
|V| = 7.2 • 10^^ GeV. The horizontal green bands correspond to the Scr region for the observed 
value for \Yb\ = (8.58 ± 0.22) x 10“^^, where we multiply for the Scr region the Icr error for 
the sake of simplicity by a factor of three. 
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will be subdominant, i.e., will lead to a reduction of the asymmetry Yi at most by a factor 
of 3. Such a reduction will still allow a generation of 1 b compatible with the observations. 
Accounting quantitatively for the wash-out effects of interest requires solving numerically the 
system of Boltzmann equations describing the evolution of the A^i, 2,3 number densities and 
of the asymmetries 11 , 2,3 in the Early Universe. Performing such a calculation is beyond the 
scope of the present work; it will be done elsewhere. 

A priori, we do not know the value of ce'''' introduced in the NLO operator in eq. (2.11). 
We will see in the next section that the low energy observables combined with Yb will give 
us information on its value. For now, as interesting cases used for illustration, we plot Yb for 
some special values of the parameters in Fig. 3. The first set/upper plot will turn out to be 
not realistic but it is still interesting because here we can see clearly, that the Majorana phases 
of the heavy right handed neutrinos are the only sources for CP violation and sufficient to 
generate Yb via leptogenesis. In this case the sign flip of Yi and Y 2 is due to the loop 
functions. The sign of sm. cj)i,i = 2,3 can be inferred from the “right sign” observation of 
Yb- The second set/lower plot is inspired by the numerical results later on. Leptogenesis is 
still successful although there c is chosen much smaller than in the first set, since we receive 
contribution from the sin 7 term, where the enhancement from the loop functions f{mi/m 2 ) 
and f{m 2 /mi) are included. Specifically, we have f{mi/m 2 ) — —2>Qf{mi/m^) for p = 7- In 
both cases Yb is dominated by Yf. The main difference between Yi and Y 2 is the efficiency 
factor; r]\Yi/'n\Y 2 — ^ for p = 7. Y 3 suffers from a strong washout in the first case and is zero 
in the second case due to cosy = 0. The NLO contribution can be regarded as an expansion 
in powers of 6 c in both cases. 

4 Phenomenology: Numerical Results 

In this section we discuss the numerical results of a parameter scan. The analytical results 
give a first impression of the general behaviour of all the observables but since there are several 
parameters involved which interplay non-trivially we made a random scan of the parameter 
space with certain assumptions to prove that our model can simultaneously fulfill all the 
constraints. The structure of this part follows the structure of the previous section. 

4.1 Masses and Mixing Angles 

For our numerical scan we follow closely the method as described in [2]. Most importantly for 
the parameters describing the quark and charged lepton sector we used the fit results given 
there. This implies that we use here tan /3 = 30 and Msusy = 1 TeV. 

In our previous model we had to scan over four real parameters (two moduli |A| and |Y|, 
two phases cj) and ^ 12 ) in the neutrino sector. In addition to these we have now scanned as 
well over the modulus c and the phase 7 . And now we have included in our scan as additional 
constraint [ 20 , 21 ] 

Yb = (8.58 ± 0.22) x 10“^^ , (4.1) 

where we multiply for the 3 cj region the la error for the sake of simplicity by a factor of three. 
For the calculation of Yb we use the formulas from section 3.2. 

Before we come to our results for the normal ordering we want to comment briefly on the 
inverted ordering. In our numerical scan we were not able to find any points in agreement 
within 3(T with all the mentioned observables. We restricted c < 0.2 and neglected points 
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Figure 4: Results of our numerical parameter scan. Blue (red) points are in agreement within 
3(j (Icr) of the low energy neutrino masses and mixings and Yb in our model. The allowed 
experimental 3a (la) regions are limited by blue (red) dashed lines. The black dashed lines 
represent the la range for the not directly measured CP phase S from the global fit [5]. 
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Figure 5: Prediction for the effective neutrino mass accessible in neutrinoless double beta 
decay experiments as a function of the lightest neutrino mass toi. The allowed experimental 
3cr (Icr) regions for the masses and mixing angles in the case of normal ordering are limited by 
blue (red) dashed lines. Blue (red) points are in agreement within Scr (Icr) of the low energy 
neutrino masses and mixings and Yg in our model. The grey region on the right side shows the 
bounds on the lightest mass from cosmology [20] and the grey region in the upper part displays 
the upper bound on the effective mass from the EXO experiment [22]. The red, straight lines 
represent the sensitivity of GERDA phase I respectively GERDA phase II [23]. 


where due to a fine-tuned cancellation the NLO corrections were artificially enhanced. Hence, 
we conclude that this ordering is still excluded like in the original model. 

The results of our scan for the masses and mixing angles is shown in Fig. 4 where the 
careful reader might note hrst that now we have as well found parameter points that are in 
agreement within Icr with all observables. That seems to be surprising since we have added 
here an additional constraint and apart from this expect rather small deviations from the 
original model. But there are two things coming together: First of all, due to the correction 
we can now allow for smaller values of 023^^® down to about 44° and furthermore we use here 
the updated results from the nu-fit collaboration [5] which allows for ^ 23 ^^® = 45° even at 
Icr. 

The second thing to note is that now the correlations between 0™^® and the phases is 
much weaker which can be explained by the fact that now we have on top another complex 
parameter in the game. But still the phases are not in arbitrary ranges but we find 

<5 G [47°, 104°] or [250°, 316°] , (4.2) 

cri E [0°,85°] or [275°,360°] , (4.3) 

cr2 E [94°,269°] , (4.4) 

For the Jarlskog invariant which determines the CP violation in neutrino oscillations we 
hnd values between ±(0.026 — 0.035). The restricted ranges for the phases imply of course 
also restrictions on the predictions for neutrinoless double beta decay, see Fig. 5. But more 
restrictive in this case is nevertheless the constraint on the mass scale where the lower bound 
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Figure 6: The relation between (j) and p according to eq. (3.14) (black line), cf. Fig. 2, and 
blue (red) points from our numerical scan in agreement within Scr (Icr) of the low energy 
neutrino masses and mixings and Yg. 


is mostly determined by the mass sum rule. We obtain for the lightest neutrino mass mi 
values between 10.5 meV to 17.6 meV. In fact, our prediction for rriee is rather precise to 
be in the narrow range from 2.3 meV to 9.2 meV. This is way below the sensitivity of any 
experiment in the near future so that any evidence for neutrinoless double beta decay would 
rule out this model. 

Related to the mass scale are as well two other observables. First of all there is the sum 
of the neutrino masses 

(0.074 - 0.089) eV , (4.5) 

which might be determined from cosmology. So far there is only an upper bound [20] 

Y'rnu <0.23 eV, (4.6) 

which is well in agreement with our prediction. The second observable is the kinematic mass 
m^ as measured in the KATRIN experiment [24] which is given as 

rrifi = miCi2Ci3 -b m2Si2Ci3 + ■ (4-7) 

Here we predict mp ~ (0.014 — 0.019) eV which is again way below the projected reach of 
mi 3 > 0.2 eV. 

4.2 Leptogenesis 

In this section we show the results of our parameter scan relevant for leptogenesis where 
we have implemented the formulas given in section 3.2 to calculate the generated baryon 
asymmetry. 

Before we actually discuss the results for the asymmetry itself we first want to note that 
the results from our analytical estimates are quite good. For instance, in Fig. 6 we show 
the relation between (p and p from eq. (3.14) and from our numerical scan. The agreement 
is striking although in the analytical estimates we have neglected for instance RGE effects 
which are nevertheless not very large in the allowed mass range. The biggest difference is in 
the allowed range for p. To avoid the resonance condition we only demanded p > 3.7 while 
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Figure 7: Results of our numerical scan for the total baryon asymmetry Yg in dependence of 
the four most relevant parameters. Blue (red) points are in agreement within 3 ct (Ict) of the 
low energy neutrino masses and mixings and Yb in our model. 


we find here p > 5.8. But here not only the ratio of the mass squared differences enter, but 
the two values of the mass squared differences independently. 

Now that we are convinced that our analytical estimates have been good we discuss the 
dependence of Yb on the four most relevant parameters as discussed in section 3.2. The 
biggest advantage of our numerical scan over the analytical estimates is that it allows us 
to use all available data on neutrino masses and mixing to constrain the allowed parameter 
space. 

We have already seen in Fig. 6 that the values of p get constrained which is again visible 
in Fig. 7. While a priori we only knew that p > 3.7 and less than about 9 we now see 
that only the range from 5.8 to 8.5 is allowed (7.4 — 8.3 at la). And since p and (p are not 
independent but related via eq. (3.14), the phase (j) gets constrained as well to the range 
[-1.2; 1.2] ([-0.5; 0.7] at la). 

The new parameters c and 7 are nevertheless more interesting than p and (j) which are 
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mostly constrained by the neutrino masses and for which we would have found similar results 
already in the previous model. In Fig. 7 we have shown the dependence of Yb on this new 
parameters. 

The first thing to note, is that c is indeed a small parameter in the range from 0.7 • 10“^ to 
4.9 • 10“^. From the model building point of view such a small value is justified. Remember 
that the leading order Yukawa coupling is a dimension three operator in the superpotential 
while the correction proportional to c is coming from a dimension seven operator. Also note 
that alone from a constraint on c could have been much larger or smaller depending of 
course on the value of 7 and the other parameters. This is different here because the mixing 
angles get corrections of order c and this implies the constraint shown here. 

Finally, note that the allowed range for 7 is only weakly constrained. Nevertheless, it is 
interesting that the sign of Yb is completely determined by 7 . This is somewhat surprising 
because in our estimates from section 3.2 the sign of depends on other parameters which 
could induce a sign flip, which can be seen for instance in the upper plot of Fig. 3. But after 
applying all experimental constraints the correlation is striking. 

Combined with the analytic analysis, we see that this correlation is a result of the fact 
that Yb is dominated by Yi, which is again dominated by the first term in ej, where cos (/>2 
is negative and f{mi/m 2 ) is positive. Neglecting the subdominant terms, we have Yb oc 
sin 7 . The analytical estimates for the efficiency factors we are using provide results with an 
estimated precision of (20-30)% compared to the full numerical results solving the Boltzmann 
equations. This is more than sufficient for the purposes of our study. The other predictions 
for the light neutrino masses and mixing parameters would only mildly change because they 
are mostly governed by the leading order values (a 30% correction to c would have only little 
impact on them). It is also worth mentioning that the complex Yukawa and the Majorana 
phases are both necessary CP-violating sources to generate a successful baryon asymmetry 
via leptogenesis while in accordance with all the low energy constraints. It is also noticeable 
that Yb would be strongly suppressed if cos (/>2 = 0. As it follows from Fig. 2 , values of 
cos 02 — 0 are excluded in the model we are considering since cos 02 can have values only 
in a narrow interval around —1, namely, (—1, —0.95) for p > 5.8, which also means that the 
Majorana phases contribute maximally to the asymmetry. In order to investigate the role of 
the Dirac phase we need a different parametrisation of the neutrino Yukawa coupling to see 
the relation explicitly, which is beyond the scope of the current work. 

5 Summary and Conclusions 

In this paper we have revised the SU(5) x A 5 golden ratio GUT flavour model from [2] 
with the aim to include as well successful leptogenesis. In the original setup this was not 
possible. As it turns out we only have to add two additional pairs of messenger fields but no 
additional symmetries or flavon fields to do this. We find that this induces a small correction 
to the neutrino Yukawa matrix, which can generate a sizeable baryon asymmetry, but as 
well implies some modifications for the predictions of the masses and mixing angles of the 
original model. In an extensive numerical scan we could show that we can simultaneously 
accommodate successfully the observed neutrino masses, mixing angles and possibly baryon 
asymmetry. And even more our setup is so constrained that we predict several correlations 
or ranges for observables yet to be measured. 

One of the most striking features of our original model - the sum rule for the neutrino 
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masses - remains valid up to a insignificant correction. From this we can again derive a 
lower bound for the lightest neutrino masses nii > 0.011 eV and rule out the neutrino mass 
spectrum with inverted ordering. This is already a very strong prediction. 

Due to the additional complex parameter and the additional constraint on Yg the allowed 
ranges for a 2 has shrunk from 70° — 290° in the original model to 94° — 269°. Whereas the 
allowed regions for ai and S remain similar compared to the original model. Namely, we find 
now ai to be in 0° — 85° or 275° — 360° and 5 to be in 47° — 105° or 250° — 316°. The strong 
correlation between 0 i 3 and the Majorana phases is now weakened due to the additional 
complex parameter we introduced in the model. It is also important to note, that we find 
here points which are in agreement within la with all neutrino observables. This is due to the 
fact that we now allow for smaller values of 6*23 but we also use here the updated fit results 
from [5] where maximal atmospheric mixing is again allowed at the Icr level. Nevertheless, a 
precise measurement of 023 which deviates significantly from maximal mixing can rule out the 
presented model. Since we limit the allowed ranges for the CP violating phases and the light 
neutrino masses we predict as well the effective Majorana mass observable in neutrinoless 
double beta decay to be in the narrow range (2.3 — 9.2) meV. This is beyond the reach of 
ongoing experiments and upcoming experiments which will begin taking data in the near 
future, but it will be certainly tested in the future. 

For the baryon asymmetry Yb we find in the approximation used to calculate it good 
agreement with the most recent data and this is done by only introducing one additional 
operator which involves one new complex parameter with a modulus c having a value in the 
range 0.7 • 10“^ to 4.9 • 10“^. The phase of this additional parameter at the 3 it level is not 
much constrained but it governs the sign of Yb- 

What we did not discuss in the present article is that some of the features of the original 
model, like the Yukawa coupling ratios yr/yb ~ ~3/2, remain valid in the modified model im¬ 
plying non-trivial constraints on the spectrum of the supersymmetric partners of the Standard 
Model particles. 

In summary we have succeeded to modify the model from [2] to include viable leptogenesis 
by only introducing a minimal correction. The model presented here is, to our knowledge, the 
first GUT A 5 golden ratio flavour model in which it is possible to have successful leptogenesis. 
All observables lie within the measured ranges and for the not yet measured quantities in the 
neutrino sector (the type of the neutrino mass spectrum, the absolute scale and the sum of 
the neutrino masses, the effective Majorana mass in neutrinoless double beta decay, the CP 
violation phases in the PMNS matrix), we make predictions. An appealing feature of the 
model is its rather small number of parameters, which makes the model very predictive and 
testable. 
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